P and T odd effects in deuteron in the Reid potential 
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The P and T odd deuteron multipoles are calculated in the Reid nucleon-nucleon potential in the 
chiral limit — > 0. The contact current generated by the 7r-meson exchange does not contribute 
to the anapole moment. The contact current generated by the vector meson exchange is negligible 
in comparison with other contributions of vector mesons. The result for the deuteron electric dipole 
moment is of great interest because of the experiment on its measurement discussed in Brookhaven. 
The deuteron photodisintegration cross section asymmetry at the threshold is also calculated. It is 
shown that its value strongly depends on the tensor forces and d-wave contribution to the deuteron 
wave function. 
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I. INTRODUCTION 

The anapole moment (AM) was introduced by 
Zel'dovich 1] as a very peculiar moment which involved 
both electromagnetic and weak interaction. A charged 
particle interaction with an AM has a contact nature. 
Actually, the interaction of electron with nuclei AM, be- 
ing of the order of qGjt , cannot be distinguished from 
radiative corrections to the weak interaction due to neu- 
tral currents. The study of Flambaum, Sushkov and 
Khriplovich 2] shows the growth of AM with the nu- 
clei size as A 2 / 3 . It means, the parity nonconservation 
(PNC) interaction of electron with heavy nuclei AM may 
become significant and provide important information 
about nuclei PNC forces. The competing contribution to 
this electron-nucleus interaction - radiative corrections 
to the neutral currents - does not have an enhancement 
in heavy nuclei. 

The deuteron AM does not have the enhancement as 
heavy nuclei and its discussion is possible for another 
reason - it has isoscalar structure only. Radiative cor- 
rections to the electron-deuteron interaction due to Z- 
exchange contains both isoscalar and isovector contribu- 
tions. The isoscalar part of interaction which is operative 
in deuteron is calculated with sufficiently good accuracy 
and is of the same order of magnitude as the inter- 
action of electron with the AM. Moreover the deuteron 
AM is defined mainly by the 7r-meson exchange and is 
singular in the 7r- meson mass. 

The study of PNC effects in deuteron has a long his- 
tory. On the one hand, the theoretical predictions are 
reliable due to small deuteron binding energy Sd ~ 2.23 
MeV. On the other hand the study of parity violation in 
the simplest nucleus - deuteron - is very important be- 
cause of existing discrepancy between experimental data 
on PNC forces in 133 Cs Q and in some other nuclei 0,01 • 

The deuteron AM was discussed in the series of pa- 
pers [2,7-11]. In [lOj the zero range approximation was 
used to calculate AM. In [ll| the deuteron AM was ob- 



tained using more realistic deuteron wave function, based 
on Argonne vig-potential. The d-wave and vector meson 
contribution in the AM were also included in the consid- 
eration. 

The same approach using effective field theory can be 
applied for P and T odd moments calculation: electric 
dipole moment (EDM) and magnetic quadrupolc mo- 
ment (MQM). The first of them is of a great interest 
because of the experiment on its measurement discussed 
in Brookhaven. The deuteron EDM and MQM have been 
considered in various models in papers [10,12,13]. 

The last part of the paper is devoted to the P odd 
asymmetry of the deuteron photodisintegration cross sec- 
tion calculation in the Reid potencial [lj]- At present, 
there is a large dispersion of the theoretical predictions 
for this value [15-30]. Most of them considered vari- 
ous models of weak PNC interaction and deuteron wave 
functions. In 28] this magnitude was calculated using 
the zero range approximation approach and the obtained 
value at the threshold was about A = 1 x 10~ 7 . The 
more complicated and reliable approaches used in papers 
|29j (with Argonne vis potential) and j^3] (with Paris 
potential) gave the following answers A = 0.253 x 10~ 7 
and A — 0.335 x 10~ 7 respectively. It seems to be very 
interesting to reveal the nature of the disagreement. 

The aim of this paper is to calculate independently P 
and T odd electromagnetic moments and the deuteron 
cross section asymmetry in chiral limit using realistic 
wave functions obtained in the soft core Reid potential. 
It is very important to reveal the factors which have the 
greatest influence on P and T odd moments: d-wave 
contribution and the deuteron wave function behavior 
at small distances. 



II. DEUTERON WAVE FUNCTIONS IN THE 
REID POTENTIAL 
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We can obtain the deuteron wave functions using the 
Reid potential Q. We represent the deuteron 3 Si — 3 D\ 



2 



state in the following form: 

<f>d = lj^(u(x) + !j=w(x)j , (1) 

where x — m^r - the dimensionless distance, Si 2 — 
3(cr 1 n)((T2n) — (er 1 or 2 ), m T = 140 MeV - charged n- 
meson mass. Normalization condition gives 

(u 2 {x) + d 2 (x)) dx = 1. 

The Shrddingcr equation for u and w components can 
be written as follows: 

u"(x) + m(E - V c {x))u(x) 

w"{x) + m(E - V c (x) + 2V t (x) + W is (x))w(x) 



m Vt{x)u{x). 



(2) 



Here V c , V t , and Vu are used for the central, tensor, and 
spin-orbit parts of the nucleon-nucleon interaction poten- 
tial, respectively: 

V(x) = V c (x) + V t {x)S 12 + Vi s (x)LS. 

The calculated functions are plotted in Figure 1. The 
obtained energy and d-wave contribution are: e = 2.23 
MeV, P d = J w 2 {x)dx = 0.065. 




x = m^r 

FIG. 1: Deuteron wave functions u,w 

It should be mentioned that these results can be ob- 
tained using either the soft core or the hard core Reid 
potential. But in the future, we will use for calculation 
of P and T effects the soft core Reid potential only. The 
hard core Reid potential is obviously inappropriate for 
very small distances description. 

III. AM NUMERICAL CALCULATIONS 

To calculate the perturbed deuteron wave function we 
will follow the way based on the direct solution of the 



Shrodinger equation. Indeed, the perturbed function can 
be written as 



Mr) 
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■jl(x) + — — i(cri + <T 2 )nu3p 1 (x) 
v8 



+i(o-i - (T 2 )nvi Pl (x) . (3) 

Functions V3 Pl and m Pl are p- waves with the total spin 1 
and respectively and the total angular momentum 1. To 
obtain these functions let us use the following weak PNC 
nucleon-nucleon potentials due to the 7r-meson exchange 

m 



V(r) 



4 7T Ul 

and p-, w-meson exchange 



99 / ^ ^ e 



(4) 



W = -9 P 



x — {{<T 1 -<T 2 ){pi-P2j P {r)}+2(l+ Xp )[<TiX(T2]Vf p (r)) 



x^((cr 1 ~ < T 2 ){p 1 -p 2 ,/ w (r)}+2(l+ Xw )[(T 1 x C r 2 ]V/ w (r)) 



-^( r i-' r 2)( cr i+ cr 2)^-{pi -P2, guhiUirygphlfpir)}. 

' (5) 
The numerical values of used parameters are listed in 

Table 1 IH. 



9p 


9u 


Xp 


X" 


h° • 10 7 


hi • 10 7 


hp ■ 10 7 


hi ■ 10 7 


hi ■ 10 7 


2.79 


8.37 


3.7 


-0.12 


-11.4 


-0.2 


-9.5 


-1.9 


-1.1 



TABLE I: Numerical values of the constants in potential 
Then the Shrodinger equation 

■^r + l l±^ + V )^ = E^ d (6) 
m r dr z mr z J 

can be split into four following equations 

u"{x) + m{E — V c {x))u{x) = 2V2 m Vtiv(x), 

w"(x)+m(E-V c {x)+2V t (x)+3V ls {x))w{x) = 2V2mV t u{x), 
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v "p 1 ( x ) ~ ~ V3 Pi ( x ) + m (E ~ V 3 Pi (x))v 3 Pi (%) 



— (F„(x)+F}(x)-F^x))- 



-1{F l p -Ft) 



-(fj, p +fx n ) j dxx I u(x) + — j^w{x) ) Vi Pl (x) 



el. 



(9) 

The numerical solution of equations gives the following 
result for the AM without nucleon contribution: 



a d 



el 



6mm T 



'14.555 + 0.048/i* - 0.132/i£ 



-(F}-F*)(u(x)-V2vj(x)), 



v'i' Pi (x) — —vi Pl (x) + m(E — Vi Pj (x))vi Pl (x) 



-0.074/1° - 0.051/i°) . 



(10) 



The nucleon AM (the pi-meson exchange only) was 
calculated in [To| 



egg 



12mm w 



6 m 
In I cr r 

7T m m„ 



(u(x) - V2w(x)) — (3 Xp F°(x) - X ,F»(x))- 



-2(3 Xp F° p (x) - XuF°(x)) — (u(x) - V2w(x)) 



'-(3 Xp F°(x) - Xu FZ(x))(u(x) + 2V2w(x)) (7) 



with used functions defined as 



Fir(x) = gg 



Attx ' 



Attx 



The general form of the AM operator is 
2tt f 

&d = Y dr[r x [r x j(r)]] + a N , (8) 

where a^v is the nucleon contribution to the deuteron 
AM. 

Let us consider the AM without nucleon contribution. 
The current operator j(x) in JSJ) can be expressed in terms 
of the perturbed wave function J2J. Then the AM is j^jj 



ad 



2n 



3m 



x / dx x(u(x) — \/2w(x))v3 Pl (x) 
Jo 



eg 



-6.19 <r 



(11) 



The last numerical result [33] was obtained considering 
vector meson exchange as well: 



12mm. 



(7.6I5 + 8.25/iJ + 2.54/i° ) Op,„. (12) 



We can see that the 7r-meson exchange contribution 
to the nucleon AM is very close to our analytical result 
which we will use further. 

To find the nucleon contribution we should average the 
sum of the nuclcons AM over the deuteron wave function: 



a N = 



u(x) + !^=w(x) ) la,, a„ ) 



( u(x) + ^w(x)^j dx = 2a p ^1 - ^P d 



Finally, the total value of the deuteron AM is 



(13) 



Sid 



6mm, 



(20.145 + 7.5/iJ - 0.132/ii 



-0.074/1° + 1.78/i° )I. 



(14) 



The contact current due to the 7r-meson exchange was 
obtained in ^| 

jc (r) = ^r(IV)^. 



2nm 

It was shown also that its contribution to the deuteron 
AM vanishes. 

The last contribution to be calculated is the vector 
meson contact current contribution. In the momentum 
representation it equals 

dW 

Jc ~ ~ ~dX • 
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But, instead of performing the calculations we can note 
that this contribution is suppressed by the factor n/m p ~ 
0.06 in comparison with other vector meson contribution 
in H14|) . It means, these values are beyond the accuracy 
of calculations. 

The total constant C2d of P-odd electron-deuteron in- 
teraction is the sum of two constants: C£ d - radiative 
corrections to the neutral current and C 2d - electromag- 
netic interaction with deuteron AM. The AM constant 



C 2d = ~ aa d 



7f 



= 0.0075 ±0.0015. 



(15) 



We have estimated accuracy on the level about 20% at 
the fixed "best value" g = 3.3 x 10" 7 [j^- 

Combining with the constant due to radiative correc- 
tions to the neutral current |3J 



C r 2d =0.014 ±0.003 



we get 



Q 



2d 



0.0215 ± 0.0035. 



(16) 



The accuracy of this calculation is high enough to give 
a hope that the experimental measurement will be able 
to provide very useful information about PNC nuclear 
forces. 



IV. ELECTRIC DIPOLE AND MAGNETIC 
QUADRUPOLE MOMENTS 

The calculation of the electric dipole moment (EDM) 
is of a great interest because of planned experiments in 
Brookehaven National Laboratory. 

Let us mention that the smallness of the vector me- 
son exchange contribution to the deuteron AM is caused 
by the smallness of the corresponding coefficients. The 
deuteron EDM and MQM matrix elements have very 
close nature and consequently, vector meson exchange 
contribution to them should also be small in comparison 
with the 7r-meson contribution. 

The deuteron EDM calculation can be performed in 
the same manner as the AM calculation was done. The 
perturbed wave function has the following form: 




+ (eri + (T 2 )nv3 Pl (x) + (<ti - tr 2 )nvi Pl (x) ) . (17) 
The deuteron EDM is 

d = (V'dlerplV'd) 



Tii r simple calculation gives the following value 

127rm w 

The deuteron magnetic quadrupole moment is 

e 



M zz = 



3mm,, 



2(/i p — pL n ) f dxxu(x)v3p 1 (x) 



2(/i p + fi n ) / dxxu(x)vip 1 (x) 



AV2 ( 

1 



dxxw(x)v3 Pl (x) 



H — -={ix p +pi n ) I dxxw(x)vi Pl (x) 
V2 Jo 



(19) 



The numerical calculation gives the following result: 

e 



M z 



1271771771^ 



(11.5sfo + 19.8 9 l 



(20) 



V. CROSS SECTION ASYMMETRY AT THE 
THRESHOLD 

Another method to observe P parity nonconservation 
in deuteron is measurement of the photodisintegration 
cross section asymmetry. Further, we will consider the 
asymmetry at the threshold of jd — * np reaction only be- 
cause it has the maximum value there. The experimental 
result of Leningrad group |34| gives for this value 



A = (1.8 ± 1.8) x 10" 



(21) 



Unfortunately, this restriction cannot give additional 
information about PNC nuclear forces. 

It is shown [2^ that this value is about 10~ 7 . But the 
value calculated is obtained in the zero range approxima- 
tion model without d-wave mixture consideration. Ac- 
cording to [18,19] these effects are significant and have a 
large influence on the final result. 

Let us consider the deuteron photodisintegration cross 
section asymmetry at the threshold. It is clear that the 
transition 3 S± -* 1 Sq is the only possible at this energy. 
We have Ml regular transition and El admixed transi- 
tion (Fig 2.). 

Only the spin nonconservation weak interaction con- 
tributes to the effect. Consequently, the deuteron wave 
functions and the wave function of the final 1 So state can 
be written as follows: 



V'd(r) 




u(x) 



S12 



w(x) + i(cr 1 - tT 2 )nvip 1 (x) 



2fi / -j^ 

/ dxx ( u(x) -I y=w(x) I V3p (x) I. (18) 

Jo V V 2 ! ' 



V^So = - if( x ) +«(o-i - (T 2 )ng{x) 



(22) 




FIG. 2: PNC transition at the threshold 



To get the value of the asymmetry we should solve the 
following system of equations: 

u"{x) + ^{E d - V c (x))u(x) = 2V2mV t w(x), 



w"(x) + ^{E d -V c + 2V t + 3V u )w(x) = 2V2mV t u(x), 



v'i Pl (x) vi Px (x) + —{E d - Vi Pl (x))vi Pl (x) 



(u(x)-V2w(x))—(3x P F°(x)- 



- Xw F° (x))- 2(3 Xp F° p (x) - x.F° (x) ) — (u(x) - V2w(x) ) 



-(3 Xp F°(x) - X.F°(x))(u(x) + 2V2w(x)), 



f"( x ) + —(E-V 1So (x))f(x) = 0, 



g"(x) -g(x) + -^(E- K Pl (x))g(x) 



^ ((2 + Xp )F° p (x) + (2 + X«)F»{x)) f(x) 



+2(F° p (x) + F°(x))f(x) - \{F° p {x) + F*)f(x). (23) 
The following symbols were used: 



F*{x)= 9uj hl 



Attx 



The asymmetry is 
A = 



4m 



3m 7r (/Lt p - j«„) J Q f(x)u(x)dx 



x ( J f(x)vi Pl (x)xdx- 



) _ 
g(x)(u(x) — V2w(x))xdx 



(24) 



The numerical value at "best values" parameters [3 
for P-odd nuclear forces is 



A = 0.16 x 10" 



(25) 



This value is much less than the asymmetry calculated 
in the zero range approximation with modified wave func- 
tions 28] . But the reasonable explanation can be given 
for this fact. The most crucial effect on the result has 
d-wave admixture to the deuteron state. Its influence 
is much stronger than the wave functions suppression at 
small distances. Actually, the same asymmetry calcu- 
lated in the Reid potential without d-wave consideration 
gives the result 

A = 0.74 x 10~ 7 . 

It means that the d-wave content in the deuteron wave 
function should be determined very precisely to calcu- 
late the asymmetry at the threshold. The Reid poten- 
tial provides P d = 6.5% for the d-wave content, whereas 
the same magnitude in Argonne v±g potential equals 
P d = 5.7%. This discrepancy can affect the result very 
strongly. 



VI. CONCLUSION AND COMPARISON 

The deuteron AM calculation based on the phe- 
nomenological nucleon-nucleon interaction was discussed 
before in papers [12,13]. The following magnitude was 
obtained 



a d — a 7r + a nucZ 



nucleoli i a p,uj 



(14.155 + 6.96g 



F°(x)=g p h° p - 



Attx 



+7.6h p - 0.14/i° - 0.2hl + 2.33/i°) I. (26) 

The 7r-meson exchange is splitted here in two terms - 
the contribution of the 7r-meson exchange between nucle- 
ons and the additive contribution of the nucleon anapole 
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moments. The 7r-meson contribution is almost the same 
in our consideration. The existed discrepancy is beyond 
the accuracy of both calculation. As to the vector meson 
exchange contributions, it differs from our result up to 
50% in some cases. But actually, this difference has the 
simple explanation - the result dependence on the d-wave 
contribution, which is different for various potentials. 
Moreover, as distinct from the 7r-meson exchange, vec- 
tor meson exchange contribution strongly depends on the 
wave function behavior at small distances (r < 0.3 fm), 
which makes the results dependent on nucleon-nucleon 
interaction model. The Reid potential, as any other 
potential description as well, fitting the scattering data 
up to 350 MeV cannot give reliable results for energy 
770 — 780 MeV needed for the vector meson exchange 
description. 

Our result for the deuteron EDM is close to that ob- 
tained in fl2| : 

d = -7^-;( 6 - 06 ffi + 2 - 37 fo + 1.05$, + 0.26s JL (27) 

The main discrepancy concerns the term with g , 
which vanishes in our calculations. It could be explained 
by the fact that paper |12| considers terms, non singu- 
lar in the 7r-meson mass. In our opinion considering this 
type terms is beyond accuracy and cannot be justified. 
The vector meson exchange contribution to the EDM is 
small, as it was expected. 

The MQM result obtained in 12] is 

M„ = -— (5.62 5o + 18.6^). (28) 



We have a reasonable agreement of our result in terms 
with gi, the discrepancy in another term is large, and 
cannot be explained because of lack of details of calcula- 
tion in 

The cross section asymmetry in — > np reaction 
was calculated in the series of papers. The latest cal- 
culations, based on the realistic deuteron wave functions 
and the "best values" constants gave the following results 

(MM)- 

A = 0.253 x 10~ 7 , 



A = 0.335 x 10~ 7 . (29) 

As we have pointed out this discrepancy can be explained 
by the strong d-wave influence on the result. Moreover, 
the above values show that the calculations of magnitude 
described by vector meson exchange as distinct from the 
7r-meson exchange cannot be reliably performed. 



I very much appreciate to V.F. Dmitriev and LB. 
Khriplovich for useful and stimulating discussions. 
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